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ABSTRACT 

The extragalactic background light (EBL) from the infrared to the ultraviolet is difficult to measure 
directly, but can be constrained with a variety of methods. EBL photons absorb 7-rays from distant 
blazars, allowing one to use blazar spectra from atmospheric Cherenkov telescopes (ACTs) to put 
upper limits on the EBL by as suming a blazar so urce spectrum. Here we apply a simple technique, 
similar to the one developed bv lSchroedteil (|2005f ). to the most recent very-high energy (VHE) 7-ray 
observations of blazars to put upper limits on the EBL energy density. This technique is independent 
of the EBL model and has well-defined errors on the constraints. Our results are consistent with EBL 
measurements and constraints but marginally inconsistent with several EBL models. 
Subject headings: gamma-rays: observations — BL Lacertae objects: general — diffuse radiation 



1. INTRODUCTION 

The EBL from the infrared through the visible and 
extending into the ultraviolet is thought to be domi- 
nated by direct starlight emission and absorption and 
re-radiation of starlight by dust. The EBL photons in- 
teract with VHE 7-rays from distant sources, such as 
blazars and gamma-ray bursts, producing e"'"-e~ pairs 
and absorbing the VHE photons. Direct measu rement of 
the EBL is difficult (see iHauser fc 0-^312001 for a re- 
view) due to contamination by foreground zodiacal and 
Galactic light. Galaxy counts may also be used to esti- 
mate the EBL, but the unknown number of unresolved 
sources results in a lower limit. Many models have been 
developed (Salamon & Stccker 1998; Malkan & Stecker' 
1998, 2001; Stecker ct al. 2006; Knciske et al. 2002, 2004; 
Primack e t allT^OOll 120 05. 200 81: IGilmore et all l2008t 
Franceschini et al.l 120081: .Razzaaue et all I2009D but the 
input physics (e.g., star formation rate, dust absorption 
and re-emission) are not constrained enough to give a 
reliable answer. 

It is possible to use VHE 7-ray observations of 
blazars from ACTs to const rain the EBL by assum- 
ing an intrinsic spec t rum JStecker fc de Jageil 119931 : 
Stanev & Franceschini' 119981 ). This was done by 
Aharonian et al. (2006a) using the hard observed spec- 
trum of lES 1101-2 32 from HESS. They ta ke a certain 
EBL shape (that of lAharonian et al.ll2003br ) that is con- 
sistent with EBL observations and lowered its normaliza- 
tion until they can fit the computed intrinsic spectrum 
with a power law softer than Tint = 1-5. This gives an 
EBL quite close to the lower limits from galaxy counts. 
This technique has, however, been criticized because it 
assumes a partic ular shape of the EBL spectrum, which is 
not well known (|Mazin fc Rauell2007f ) . This criticism ap- 
plies to severa l other constraints from b laz ar o bserv ations 
as we ll (e.g., lAharonian et al.l l2007d^ . "Mazin fc Raud 
(|2007f ) have addressed this problem by developing a so- 
phisticated technique which scans over a large grid of 
possible EBL shapes. They de-absorb blazar VHE 7- 

Electronic address: [Justin. finke®nrl.navy.iml| 



ray spectra with every possible shape, and throw out the 
shapes which, when fit with a power-law, broken power- 
law, or power-law with an exponential cutoff, give Tint 
greater than a certain value (they get two constraints, 
from Tint > 1.5 and > 0.75). The number of remaining 
allowed shapes give an upper lim it on the EBL intensity . 
A similar technique was used bv lKrennrich et all (|2008f ) 
to rule out EBL models which ar e inconsistent with 
the r ecent lower limit at 3.6 /zm (jLevenson fc Wrightl 
2008f). A relative ly simple technique was developed by 



Schroedtej (I200I to give upper limit s on the EBL which 



is simpler than the techniqu e of ,Mazin fc Raui ()2007f ) 
and has well-defined errors. ISchroedteil (|2005f ) applied 
his technique to the 6 blazars seen with ACTs at that 
time to get upper limits on the EBL energy density. To 
date, VHE 7-rays have been seen from 21 blazars, in- 
cluding 1 flat spectrum radio quasar (FSRQ) and 20 BL 
Lac object s . In this paper, we apply the technique of 
iSchroedtej (I200I to an up-to-date sample of TeV blazar 
spectra. 

Note that this technique, and other techniques to con- 
strain the EBL with 7-ray observations from blazars 
do make assumptions about the intrinsic spectrum of 
blazars, and is thus not a limit in the conventional sense. 
That is, understanding of blazars' intrinsic spectra could 
improve in the future, which would change our conclu- 
sions. If blazars were found to typically have harder 
spectra than we assume here, our upper limits would 
be weaker, whereas if they were found to have typically 
softer spectra, they would be stronger. 

In § 2 we describe how one can put upper limits on the 
77 absorption optical depth, r^^ (e) , and the EBL energy 
density, euEB hi^^ z), usiiig a me thod similar to the one 
developed by ISchroedteil ()2005l ) . In § 3 we apply this 
technique to the most recent ACT spectra from blazars 
and discuss our results. 



2. CONSTRAINING THE EBL ENERGY DENSITY FROM 
GAMMA-RAY OBSERVATIONS 

The observed vF^, spectrum of a distant source, fobs, 
at redshift z is related to its unknown intrinsic spectrum, 
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fint by 



fobs{e) 



/mi(e) 7 



(1) 



where e is the dimensionless observed 7-ray photon en- 
ergy. The greater the absorption, the harder the intrinsic 
spectrum must be to produce a given observed spectrum. 
The hardest spectrum observed in the TeV range from a 
blazar near enough not to be strongly affected by EBL at- 
tenuation is r fti 1.9 (where /(e) = .f (e m»n)(£/emm)~^'^^) 
from Mrk 501 ('A haronian et al.lll999l ). Two blazars ob- 
served with EGRET have T « 1.5, the FSRQ 0847- 
120 and the TeV-observed high-peaked BL Lac Mrk 
501 (|Nandikotkur et all [20071) . The hardest blazar with 
the Large Area Telescope (LAT) in the Fermi Gamma- 
Ray Space Telescope 3-month bright AGN source list is 
r ^ 1.4 from the B L Lac MS 1050.7+4946 at z = 0.140 
(|Abdo et al.l I2009D . This blazar has not been seen in 
the TeV range, however, and may be a different class of 
blazar. Just over half of the TeV blazars have been seen 
with Fermi and all of the TeV High-Peaked BL Lacs have 
Fermi spectral indices F between 1.70 and 1.85. From 
a theoretical standpoint, shock- accelerated electrons are 
unlikely to produce 7-rays with spectral indices harder 
than F = 1.5 from Compton scattering. We will there- 
fore conservatively take F^t > F™™ = 1.5 as the con- 
ventional limit on the spectral index. Several theoretical 

fjssibiliti es have b een proposed to create harder spectra. 
atarzvnski et al.l ([2006.) have suggested that electrons 
create synchrotron and Compton emission at a signif- 
icant distance from the acceleration region, leading to 
a high cutoff in the lower portion of the electron spec- 
trum that generates the TeV 7-rays. This could result 
in harder observed spectra. IStecker et al.l ()2007f ) have 
performed simulations which generate harder electron 
spectra at the highest energies in relativistic shocks, and 
hence harder observed TeV 7-ray spectra. Internal 77 
absorption could possibly produce harder spectra in the 
TeV range (jAharonian et al.|[2008bl ). although this would 
probably apply more to distant FSRQs with significant 
scattered radiation fields (i.e., broad line regions) which 
we do not consider in this paper. The decay of pions from 
a hadronic source could also produce an extremely hard 
TeV component, indepen dent of the lower e nergy elec- 
tron synchrotron emission (|Mucke et al.ll2003D . Compton 
scattering of the cosmic microwave background radiation 
in t he extendedjets of bla zars could produce harder spec- 
tra (jBottcher et al.l [20081 ). These are possibilities that 
must be considered. Furthermore, creating a model con- 
sistent with the recent gal a xy count result a t 3.6 /x m by 
iLevenson fc Wrightl (|2008f) . iKrennrich et al.l (|2008f) find 
that spectral indices harder than Tint ~ 1-5 seem neces- 
sary. We will, therefore, take Tint < — l-^ as the 
extreme limit. 

By assuming the intrinsic spectrum is limited by F™™ 
we can put an upper limit on r^^(e). 
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fint (^min ) 



_ J obs y^min J 

(F„fc,-F™r)ln(e/e„, 
The standard error of rj^^^^^ is given by 

CT(/o6s(e))//ofcs(e) . 



+ 
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(3) 



We use the 7-ray flux at Emin = ^mm'mec'^ TeV, the low- 
est energy bin of the observation, to normalize the ob- 
served spectra. We then use the EBL model which gives 
the greatest T^-,(e min) (which is almo st always the fast 
evolution model of lStecker et al.l ()2006l )) to de- absorb the 
spectra here, i.e., frnt{f-m.%n) = /ahs(emm)e'"-^^^'^'"'"\ and 
use this to normalize our intrinsic maximum spectra. 

Once T"!y°'^{e,z) is found, one can calculate an upper 
limit on the local EBL. This is done using several ap- 
proximations, as follows. We begin with the 77 opacity 
of the universe as a function of the comoving EBL energy 
density, which is given by 
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{i + z'Y 



dtt 
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de' 



'■EBL 



(6';z') 
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where, in a flat ACDM universe, 

dt^, 
d? 



(6'ei(l + z')) (4) 



(5) 



Hoil + Z')^n,n{l + Z'Y+^A 



''EBL 



(e'; z') is the comoving EBL energy density per unit 



comoving dimensionless energy, e' = e(l + z), 

<P('So) -j ^ In Wo - Po In Wo - z ^ 

i — Pq ^ — Po 

+2(3q + 41nwQ ln(l -I- wq) - 4L{wo) , 
/32 = 1 - 1/so, wo = {l+ /3o)/(l - /3o), and 



L{wo) 



dw w ^ ln(l + w) 



(6) 



(7) 



(jCould fc Schredeilll967l : iBrown et al.lfl973l) . 

If we assume the blazar is at low redshift, then 
\dt*/dz\ « Hq^, where we use Hq = 70 km Mpc~^ 
for the Hubble constant. We also use the approxima- 
tion that the EBL does not change signiflcantly at low 
redshifts, and that all of the EBL photons are at the en- 
ergy where the pair-production cross-section is largest, 
e'^ = 2/(ei(l -\- z')) w 2/ei (i.e., we assume the EBL 
is monochromatic). The latter approximation overesti- 
mates the actual EBL energy density, since in reality the 
absorption will be spread out over a range of EBL pho- 
ton energies, so that we are left with an upper limit on 
the EBL energy density, e' u'^^f^ {e; z) . Using the Dirac 
delta-function these approximations are written as 



,'gr«;^«0) Sie'-e'J .(8) 



"£;sL(e';^ « 0) w e'^u'j^BL 

Eqn. ([5]) allows us to perform the integrals in eqn. ^ 
so that 



'(euz^O) 



I I max ( J . ~ , 
, ^*"-EBL z 
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or 



<««;^«o)^ 



3c(TTZ0(2)ei 
xr7^'^-(2/6:,z«0) 



(10) 
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where 0(2) w 1.787. The error in this EBL hmit is given 

by 



/ frnax\ 



^Kr) (11) 



3ccrT(/'(2)zei 

where (j{tJ^^), the error in t"1^^ , is given by eqn. ([3]). 
For z w 0, the comoving EBL energy density will be 
equal to the observed energy density, i.e., 

e*WBBL(e*,-z w 0) = f-WEBL{<^'*iZ w 0) . (12) 
For higher redshifts, 

e*UEBL{e*; z) = {1 + z)~^('_^ u'EBLi<^U^) ■ (13) 

The EBL energy density can be converted to intensity in 
units of, e.g., nW sr~^ by 



eleiz) = — t^UEBL{^*\z) 



(14) 



By inspecting eqns. ^ and PU)) we see that 
eu^1fj^{e; z) will be minimized, thus giving the strongest 
constraint, if the dimcnsionless parameter 



TeV \ Tobs - 



In 



En 



(15) 



is minimized, where Emax = "meC^^max is the energy of 
the highest energy photon bin observed from a source. 
This parameter should be seen as a general guide. The 
highest energy photon bin is often not well observed, for 
example, in H 1426 +428 the highest energy bin is less 
than a 2a detection (jAharonian et al.ll2003bh . 

3. RESULTS AND DISCUSSION 

The blazar sample used to constrain the EBL is seen 
in Table [TJ This is not a complete sample of TeV blazar 
spectra^ , as we have selected blazars that are nearby and 
should give the strongest constraints, i.e., that have the 
smallest ^. The EBL upper limits from the 8 blazar 
spectra with the lowest ^ are plotted in Fig. [1] 

For each blazar, i, one gets a list of eu^lf^ ^(e) at each 
EBL energy bin j, which corresponds to an observed 
7-ray energy bin. For N such blazars, we choose the 
min{eu^^^ j (e)} to give the strongest constraint on the 
jth energy bin. These constraints, along with other EBL 
measurements and constraints, are plotted in Fig. [21 

In section § [5] we have made several assumptions to 
get the upper limit, namely, that the EBL does not 
evolve at low redshift, and the absorption occurs en- 
tirely from an EBL photon of energy equal to the en- 
ergy where the absorption cross-section is maximized. 
There is no direct way to measure the EBL at higher red- 
shifts, however the star forma tion rate increases between 
z = 0.0 and z = 0.2 (e.g., [Hopkins fc Beac om 2006') 
as d oes the luminosity density at various wave lengths 
Ce.g.. lLe Floc'h et al.l[2005l : IBabbedee et al.ll2006D . Thus 
it seems logical that the EBL energy density would be 
greater at higher redshifts. A greater EBL energy density 
in the past would result in a greater absorption, which 
would give weaker upper limits than what the actual lim- 
its would be if evolution was taken into account. Thus, 

^ For complet e, continually upda ted catalogs see the TeV- 
Cat ( http : //tevcat . uchlcaRO . edu/ ) or R. Wagner's website 
((httpjV/www . mppmu . mpg . de/ ~ rwagner / sour ce s/ ) . 



if we did take evolution into account, our limits would 
be stronger. But the largest source of error is likely to 
be the monochromatic assumption. If the EBL energy 
density is significantly greater at higher energies (shorter 
wavelengths) than where the cross section is maximized 
(eie' = 2), the result is a considerably weaker upper limit. 
In all of the EBL models presented in Fig. [H the EBL is 
considerably higher at > 0.2 eV (< 5 /xm) than immedi- 
ately below this energy. Thus, the EBL energy density 
at energies less than 0.2 eV is likely to be considerably 
lower than our upper limits, although to actually quan- 
tify this error one would need to assume an EBL model. 
Since our conservative assumptions all result in upper 
limits erring on the high side, our constraints are quite 
strong. 

These EBL upper limits ma y be c ompared to the 
constraints by [Mazin fc Raud (|2007f ). Their tech- 
nique gives similar results and is in some ways com- 
plementary to our results. We use more constrain- 
ing blazars than iMazin fc Raud (|2007f ) , particularly the 
blazar lES 0229-1-200, which gives a stronger constraint 
at ^ 0.04 eV (^ 30 /xm). The repo rted spectrum of 
lE S 1101-232 (lAharonian et al.l l2006al). which was used 
by [Mazin fc Raud (|2007[ ). ha s been re-analyzed wit h 
improved energy calibration ([Aharonian "eTal] [2007bl ). 
which we use in stead. This is one reason our constraints 
are weaker than [Mazin fc Ra uc (2007) in the 0.2-0.5 eV 
(3-7 /im) range, azin fc Raue, (20070 also find that 
the size of the EBL grid to which they fit splines is a 
cause of systematic error. By using a finer grid, they 
find they can get constraints 20-30% higher. This could 
be another source of discrepancy between their results 
and ours, since our results do not make any assumptions 
as to the coarseness of the possible shape of the EBL 
spectrum. Additionally, our monochromatic assumption 
could lead t o a Jiigher upper limit, as discussed above. 
[Mazin fc Raud (|2007h estimate their errors in the optical 
and near-IR to be about 30%. Our results are within 
these errors at 0.2-0.5 eV. Our results are lower in the 
optical (< 10 /im or > 1 eV) due to the assumption of 
the highest possible EBL energy density. They choose 
an arbitrary level larger than we do, since we choose 
the largest p o ssible model, the fast evolution model of 
[Stecker et al.[ ([2006f ). The choice in this region does not 
act as a constraint, since our results will naturally be 
above the highest EBL model. 

Our results are consistent with the recent highly 
constraining EBL lower limit at 0.34 eV (3.6 fj,m) 
([Levenson fc Wright! [2008f) from number counts with 
the Spitzer Space Telescope (see [Krennrich et al.[ [2008[ . 
for a different int e rpreta tion). The lower limit of 
[Levenson fc Wright[ ([20081 ) combined with our results 
and other data means that the EBL is quite well con- 
strained at 3.6 /im. T he only mode l abov e this limit 
is the best fit model of [Kneiske et al.[ ([2004[ ). However, 
their model falls just above our conventional upper limit 
(r pr = 1.5) at » 0. 7 eV (p^ 1.7 /xm). The models 
of lCilmore et al.[ ([2008[ ). [Franceschini et all (12008) . and 
the fast evolution model of [Stecker et al.[ ([2006^" " also fall 
above our conventional limi t at ~ 0.04 eV (~ 3 /im), 
while the baseline model of [Stecker et al.l ([2006f ) is just 
below our upper limit. None of the models are above the 
extreme limit (F™™ — 1), which is not surprising, con- 
sidering this weaker limitation was designed to avoid con- 
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Fig. 1. — Upper limits for the conventional (F™™ = 1.5; black filled inverted triangles) and extreme (F^j" = 1.0; red empty inverted 
triangles) limits on the EBL. A lso plotted are several EBL models: the best fit m odel from 'Knciske et al.l l|2004l . solid green curve), 
the fiducial model from [Gilmore ct al. (20 081, dotted blue curve), the model of IFran ccschini ct al. (2008, short dashed violet curve), the 
baseline and fast evol ution models of Stecker et al.l (120061 . lower and upper long dashed brown curves, respectively), and model B from 
I Razzaque et al.l (120091 . dot-dashed orange curve). 

flicts between 7-ray observations and EBL models. Fu- 
ture observations of blazars with ACTs as well as the re- 
cently launched Fermi Gamma- Ray Space Telescope will 
provide further upper limits on EBL models. 

We have updated the technique of Schrocdteil (j2005l ) 
to a recent sample of VHE 7-ray blazar spectra to obtain 
upper limits on the EBL energy density. This technique 
does not make assumptions about the shape of the EBL 
or fit a certain model to a de-absorbed spectrum, thus 
avoiding problems with assuming an EBL spectral shape 
in other techniques for using VHE 7-ray observations to 
constrain the EBL. We obtain upper limits weaker and 
more conservative than limits from de-absorbing spectra 
with an assumed EBL shape and fitting the results with 



a power-law (e.g.. lAharonian et a l. 2006a) and in agree- 
ment with a similar technique bv ,Mazin fc Rau3 (|2007| ). 
Our conservative assumptions imply that it is likely the 
actual EBL energy density is considerably lower than our 
upper limits, especially for meC^e ^ 0.2 eV. 
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